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Abstract 

The frequency noise power spectral density of a free-running quan¬ 
tum cascade laser frequency comb is investigated. A plateau is observed 
at high frequencies, attributed to the quantum noise limit set by the 
Schawlow-Townes formula for the total laser power on all comb lines. In 
our experiment, a linewidth of 292 Hz is measured for a total power of 
25 mW. This result proves that the four-wave mixing process, responsible 
for the comb operation, effectively correlates the quantum noise of the 
individual comb lines. 


Introduction 

In recent years, optical frequency combs (OFCs) have become fundamental tools 
for near infrared (NIR) spectroscopy and metrology [U |2]. Thanks to their 
wide spectral coverage, high coherence and absolute traceability, they are used 
for high resolution and precision atomic and molecular spectroscopy in this 
spectral region [3]. Since the fundamental ro-vibrational transitions of simple 
molecules fall in the mid infrared (MIR), it is of particular interest to have 
OFCs operating in this spectral region. Until now, a well-established approach 
to satisfy this need consisted of directly transferring NIR OFCs emission to the 
MIR region through nonlinear processes, e.g. difference frequency generation 
using intense fiber-based NIR OFCs HSIS] or optical parametric oscillators 
[3 ig. This approach guarantees good spectral coverage and coherence, but 
requires complex and delicate experimental set-ups. Comb generation has also 
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been achieved by parametric oscillation in high-Q micro resonators [S], with a 
spectral coverage recently extended to the MIR region [13 E]. 

Quantum cascade lasers (QCLs) [12] are current-driven semiconductor lasers 
based on intersubband transitions in quantum wells, emitting MIR or terahertz 
radiation. In devices designed with low group velocity dispersion, it has been 
shown that comb operation can be achieved [HEll thanks to the four-wave 
mixing process taking place in the gain medium itself m- For MIR-operating 
devices, the upper state lifetime, inherent to the intersubband transition of the 
active region, is very short (sub-picosecond range). This is responsible both 
for the broadband nature of the four-wave mixing process that enhances the 
mode-locking, but also for a tendency to operate with a mostly constant output 
power. For these reasons, the phase relation between the modes is similar to 
that of frequency-modulated lasers m, as predicted by theory [lads, and no 
pulses are emitted. 

Quantum cascade laser frequency combs (QCL-combs) were initially char¬ 
acterized by measuring the autocorrellation of the intermode beat note at the 
cavity round-trip frequency (7.5 GHz), performing a so-called beat note spec¬ 
troscopy [T3]. A more sensitive technique is provided by comparing two QCL- 
combs in a heterodyne beat experiment. Recent experiments on QCLs in a 
dual-comb spectroscopy set-up demonstrated a mode equidistance as good as 
7.5 X 10“^® relative to the carrier optical frequency a value close to those 
measured for micro-resonator-based combs [Sj- 

In this manuscript we investigate the frequency noise of such combs. This 
characterization is essential both for spectroscopy applications as well as for a 
better understanding of the fundamental properties of these devices. The gen¬ 
eration of the comb of frequencies is interpreted within the framework of super¬ 
modes. A high-finesse optical cavity has been used as a multimode frequency- 
to-amplitude converter to retrieve the QCL-comb intrinsic linewidth. A com¬ 
parison between the linewidth obtained in the comb regime operation and in the 
single-mode operation is also given, demonstrating that the four-wave mixing 
process effectively correlates the quantum noise of the comb lines. 

1 Frequency noise 

What distinguishes a frequency comb from a simple array of perfectly equally 
spaced single-frequency optical sources is the correlation of the frequency noise. 
While the heterodyne beat of two independent single-frequency laser sources 
always yields a linewidth wider than that of the individual lasers, this is not the 
case if the two single frequencies are extracted from a frequency comb source. 
These considerations are equally true for technical as well as for quantum noise. 
The intrinsic linewidth of a laser is given by the Schawlow-Townes formula 
m and can be interpreted as the ratio of the number of photons emitted in 
the cavity by spontaneous emission over the number of the ones emitted by 
stimulated emission. As compared to a single-frequency device, the only effect 
of comb operation is the redistribution of the stimulated photons into equally 
spaced modes. For this reason, the intrinsic linewdith of a single comb line 
is expected to be unchanged and can be expressed by the Schawlow-Townes 
formula considering the total optical power of all comb lines. 

As in micro-resonators-based combs, QCL-combs are generated through a 
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four-wave mixing process. For this reason, starting from a semiclassical ap¬ 
proach for QCL-combs HU ST], it makes sense to compare the quantities with 
the quantum formalism developed for micro-resonators-based combs |19) to re¬ 
trieve the Langevin equation for the photon annihilation operator related to the 
n-th QCL-comb line: 
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with K,n,k,i,m Spatial superposition integral among the modes involved in the 
four-wave mixing process, photon lifetime in the laser cavity, 722 scattering 
rate out of the excited laser state, 712 loss of coherence of the laser transition, 
/rep comb line spacing (cavity round trip, without dispersion), peak gain, 6n 
difference between the frequency of the n-th mode of the ideal laser cavity /„ and 
the frequency of the n-th mode of the laser cavity with dispersion. Vn are the 
vacuum Langevin noise operators related to the optical loss processes (waveguide 
and mirrors ) [iniiis], characterized by the following statistical properties: 


{Vn{t)) = 0 

{V^{t)Vr,{t'))=5{t-t') 
{Vu{t)Vn{t')) = 0 


The second term in equation is due to the four-wave mixing and it is responsi¬ 
ble for the coupling among all the laser modes. Solving equationj^is beyond the 
scope of this paper. However, we note that in comb operation the average rel¬ 
ative phases of the modes are fixed (within fluctuations). Therefore, through a 
unitary transformation, it is possible to select a new basis for the cavity modes, 
the supermodes basis, such that one of this supermodes corresponds to that 
selected by the comb operation m- The new annihilation operators are given 

by 

h, =Y,U^ dn (3) 
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where is the element of a unitary matrix such that [52]. In this 

way, the equation is reduced to that of a single-mode laser (with only one mode 
excited). In particular, the Langevin operators for the new modes 

(4) 

n 

will have the same correlation properties (equations because of the unitary 
nature of the transformation. The resulting frequency noise is expected to be 
the same of a single-mode laser. 


2 Experiment 


The laser used for these experiments is a QCL-comb based on an InGaAs/InAlAs 
broadband design with multiple active regions (multistack), previously reported 
in m- It operates in continuous wave at room temperature emitting several 
milliwatts of power at 7.10 /rm. The device length is 6 mm, corresponding to 
/rep ~ 7.5 GHz. Two main operation regimes are observed in this device. Just 
above threshold, the device emits single-mode radiation. A comb regime is then 
observed for a significant part of the device working range (see figure . This 
very convenient property enabled us to study and compare the linewidth in both 
regimes (single-mode and comb) in the same device. In order to investigate the 
frequency noise power spectral density (FNPSD), a high-finesse optical cavity 
{Fabry-Perot, T « 6000) was used to resolve the laser spectrum and to detect the 
frequency fluctuations of the laser, acting as frequency-to-amplitude (FA) con¬ 
verter (figure [^. In the set-up, an optical isolator (transmission T = —2.9 dB, 
extinction E = —33 dB) has been used to avoid the instabilities induced on 
the laser by the back-reflection from the input mirror of the cavity. To collect 
the signal transmitted by the cavity, a high-sensitivity nitrogen-cooled HgGdTe 
(MGT) detector {BW = 0—10 MHz) was used. A 12-bit vertical resolution, 
1 GHz analog bandwidth, 2.5 GS/s sampling rate oscilloscope was used to col¬ 
lect the signal and to compute its Fourier transform. The distance between the 
two mirrors was chosen in order to set the free spectral range {FSR) of the cav¬ 
ity close to the comb line spacing /rep. In order to resolve the laser spectrum, 
a Vernier ratio Vr = FSR/frep slightly different from one was chosen and a 
piezoelectric actuator was used to scan the cavity length over one FSR [23]. A 
schematic representation is depicted in figure |^. 

To utilize the cavity as a FA converter, we acted again on the piezoelectric 
actuator and on the temperature controller of the laser to set Vr = 1.0 and 
to let the comb offset frequency /□ be equal to the one of the optical cavity. 
In this way, the comb lines and the optical cavity resonances were perfectly 
matched (see figure §>)■ As a consequence, in these conditions and only in 
these conditions of temperature and driving current of the laser, all the comb 
lines were transmitted by the cavity. The cavity can thus be used as a mul¬ 
timode frequency-to-amplitude converter to collect the frequency fluctuations 


of all the modes at the same time |6] (see Supplement 4.1 for a demonstra¬ 
tion). Since Vr = 1.0, an accurate value of the optical cavity free spectral 
range FSR can be obtained by measuring /^gp. The latter was measured by 
acquiring the intermode beat note on the laser current with a spectrum ana¬ 
lyzer [23]. An accurate value of the cavity FSR is needed for the calibration of 
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Figure 1: Top: Experimental set-up used to measure the FNPSD of the laser. 
The main optical components include the laser (a multistack luGaAs/InAlAs 
QCL), the optical isolator, the high-hnesse optical cavity and the high-sensitivity 
MCT detector. The signal is processed by a high-sampling-rate oscilloscope. 
Bottom: Power-versus-current curve of the QCL at fixed temperature. Two 
operation regimes are observed in this device, a single-mode regime and a comb 
regime. 
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Figure 2: Schematic of the optical cavity and of the comb spectra with Vr ^ 1.0 
(comb regime - one mode) and Vr = 1.0 (comb regime - all modes) (a) and b) 
respectively). 

c) Cavity transmissions acquired in the three conditions: single-mode QCL, 
QCL in comb regime with only one mode in resonance with the cavity, and 
QCL in comb regime with all the modes in resonance with the cavity. These 
acquisitions are obtained by scanning the cavity length. The cavity detuning is 
the variation of the resonance frequency (FSR) with the length. These peaks 
are used for the calibration of the FA converter. 
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the FA converter. Under optimal conditions of current and temperature for the 
matching, the laser emits a power P = 25 mW. Thanks to the high finesse, it is 
also possible to collect the frequency fluctuations of an individual comb line by 
slightly varying the FSR. A spectrum retrieved with the laser in single-mode 
operating conditions {P = 15 mW) was also acquired. The FNPSD measured 
on the single-mode regime as well as on the comb regime are reported in fig¬ 
ure after correction for the finite bandwidth of the FA converter (response 
frequency cutoff of the high-finesse cavity related to the photon lifetime). We 
observe that the frequency noise on the comb regime is identical to the frequency 
noise on the single-mode regime. Moreover, the frequency noise of an individual 
comb line is also equivalent to the one acquired on all the comb lines together. 
By integrating the FNPSD using Elliott’s formula [IS], the full width at high 
maximum (FWHM) of a laser mode can be retrieved. In our case, we obtain 
a FWHM of about 600 kHz in a 1-s timescale, which is consistent with the 
linewidth shown by distributed-feedback (DFB) QCLs [MlEllEH]- Moreover, 
the contributions of the driver noise and the laser intensity noise as well as the 
detection noise floor are reported. Taking into account the detection noise floor 
shape, the spectra are reliable up to 3 MHz. Around 1 MHz, a flattening can 
be observed. Figure shows a portion of the same FNPSD spectrum (from 
100 kHz to 3 MHz). This flattening, characteristic of a white frequency noise, 
corresponds to the intrinsic quantum noise level Dsi, due to the spontaneous 
emission, the so called Schawlow-Townes limit m- 


3 Discussion 

It is interesting to compare this level Dsi, to the one expected for a single-mode 
device with the same characteristics 1^: 


hv atotc^ 
P Airnl 


m ' t-sp 
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Taking i/ = 42.2 THz as central frequency, P = 25 mW as total power emitted 
by the laser, = 2.2 cm“^ as mirror losses, atot = 7.2 cm“^ as total losses (the 
relatively high waveguide losses are due to the residual cross-absorption given by 
the multistack structure), rig = 3.4, rigp = 2 as spontaneous emission factor and 
(Og) = 0.0023 as squared Henry linewidth enhancement factor averaged over 
the laser spectrum (see Supplement 4.21, we obtain Siy = 230 Hz, value which is 
consistent with the one obtained from the spectrum Si/ = ttDsi, = (292 ±78) Hz 
(figure]^). This justifies the theoretical framework introduced in section]^ 
More importantly, the measurement of the FNPSD in comb regime shows 
that the quantum fluctuations of the different modes are correlated. In fact, 
we observe that the FNPSD - in particular the portion limited by the quantum 
noise - is identical when measured with one comb line and with all comb lines 
simultaneously. This quantum limit - value which is given by the Schawlow- 
Townes expression - would be at least a factor of 6 larger than the one showed 
in figure [^, assuming that the quantum fluctuations of each comb mode are 
uncorrelated. This factor is outside our uncertainty. 
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Figure 3: a) FNPSD of the QCL-comb taken in three different conditions. The 
spectra are compensated for the FA converter cutoff. The technical contribu¬ 
tions to the noise are also reported: taking into account the detection noise floor 
shape, the spectra are reliable up to 3 MHz; the two contributions given one by 
the current driver and the other one related to the intensity noise are negligible, 
b) Zoom of the flattening portion of the spectrum around 1 MHz, corresponding 
to the Schawlow-Townes limit. 





















































Conclusions 


With this work we have demonstrated that in quantum cascade laser frequency 
combs the four-wave mixing process at the origin of the comb operation regime 
is also responsible for the correlation between the frequency fluctuations be¬ 
tween the modes until the quantum limit. As a result, the linewidth is shown 
to be limited by the Schawlow-Townes formula, as it is for single-mode lasers of 
the same total power. As a consequence, instruments using the spectral multi¬ 
plexing of dual-combs or multi-heterodyne spectrometers hold an inherent noise 
advantage compared to similar systems using arrays of single-mode lasers. Fi¬ 
nally, the same technique used to retrieve the frequency noise power spectral 
density could be used to implement an active stabilization for locking this combs 
to high-finesse ultra-stable optical cavities. 
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4 Intrinsic linewidth of quantum cascade laser 
frequency combs: supplementary material 

4.1 Demonstration of the multimode frequency-to-amplitude 
conversion 

In this section we will show how satisfying the right conditions it is possible to 
use an optical cavity as multimode frequency-to-amplitude converter to retrieve 
the frequency noise of an optical frequency comb. The intensity of the z-th comb 
line transmitted by the optical cavity is given by 

li = ( 6 ) 

where lio is the intensity of the Tth comb line, Vij = Vj — Vi, Vi is the center 
frequency of the i-th comb line, Vj is the center frequency of the j-th cavity 
transmission and T[vij) is the convolution between the i-th comb line shape 
and the j-th cavity transmission profile. The total intensity transmitted by the 
cavity is given by 

Itot = ^ h = ^ [hoTivij)] 

i i 

If 

Vij = Av = const. Vi,j 
we are allowed to bring T{v) out of the sum, yielding 

/tot = T{Av) Y, ho = T{Av)Io tot 

i 

This proves that in this conditions the transmission and consequently the frequency- 
to-amplitude conversion follows the same rules of the single-mode case. Actually, 
the condition expressed by equation which corresponds to impose that both 
the comb lines and the cavity resonances have the same dispersion in frequency, 
is not so strict. It is in fact sufficient that 

« ^FWHM [r(Aj.)] (10) 

for all the transmitted peaks. In other words it is sufficient that the accumulated 
dispersion over the whole comb spectrum is small compared to the width of the 
cavity resonances. 


(7) 

( 8 ) 
(9) 


4.2 Henry linewidth enhancement factor estimation 

In order to compute the Schawlow-Townes limit 

hv Q-totC /ii2\ 

-kae) (H) 

a proper estimation of the Henry linewidth enhancement factor is needed [29j . 
As reported in |30| . assuming a Lorentian lineshape for the intersubband tran¬ 
sition, the Henry factor for a QCL can be expressed as 


ae{v) 


V32 - V 
732 


( 12 ) 
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where vz 2 and 732 are the center frequency and the half with at high maximum 
of the transition respectively, while v is the frequency of the radiation emitted 
by the laser. Concerning we have to consider that the emitted radiation is a 
comb of frequencies, so Og has to be computed as average over the whole emitted 
spectrum 

(Qie)laser spectrum — “7 ^ / (^32 ^ ) dr' (Id) 

Al/ 732 dlas. sp. 

where Av is the width of the comb spectrum. Assuming that the center of the 
emitted spectrum corresponds to 1^32 we obtain directly 


/laser spectrum 



(14) 


Finally, knowing from electroluninescence measurements that 732 = 150 cm“^ 
and that Az/ = 25 cm“^ in the actual operating conditions of the laser (P = 
25 mW) we have 

(^e)laser spectrum — 0.0023 (15) 


value that can be neglected in equation 11 In 


shape of the gain curve on the top is even flatter 
the given value for (a^) has to be considered as 


case of a multistack device the 
than a Lorentian, consequently 
an upper limit. 
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